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Abstract. In this article we outline the structure of a general relativistic astro- 
metric model which has been developed to deduce the position and proper motion 
of stars from 1-microarcsecond optical observations made by an astrometric satellite 
orbiting around the Sun. The basic assumption of our model is that the Solar System 
is the only source of gravity, hence we show how we modeled the satellite observations 
in a many-body perturbative approach limiting ourselves to the order of accuracy 
of (v/c) 2 . The microarcsecond observing scenario outlined is that for the GAIA 
astrometric mission. 
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' Abbreviations: GR - General Theory of Relativity; CM - Center of Mass 

O ' 1. Introduction 

c3 Aim of fundamental astrometry is to measure stellar positions, pa- 

rallaxes and proper motion directly from optical observations. Under- 
standing how electromagnetic waves propagate in a time dependent 
rN \ gravitational field is essential to answer basic questions in astrophysics. 

In fact, most of the physical information about the surrounding universe 
(see Kopeikin et al., 1999, Kopeikin, 2001 and references therein) is 
carried by light signals. Nowadays the scientific objectives of astrometry 
and astrophysics can be combined thanks to the precision achievable 
by the new generation of space astrometry missions like GAIA of the 
European Space Agency (ESA). GAIA will be launched not later then 
2012 (with a possible window in 2010) and is expected to make as- 
trometric measurements with a precision of a few /iarcsec (Perryman 
et al., 2001). At this level of accuracy one should consider contributions 
to the light deflection not only by the mass of the Sun and of the 
other planets, but also by their gravitational quadrupole moments as 
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Table I. Detectable first order post-Newtonian relativistic 
effects in Solar System at 1 /iarcsec level. Here ^XpN is the 
contribution due to the spherically symmetric component of 
the gravitational field of the corresponding body, J2 is that 
of the quadrupolar moment, and L that of the rotational 
motion. The values whitin the parenthesis represents the 
maximum angular distance between the perturbing body and 
photon at which the effect still attains 1 /uarcsec; where two 
values are reported, they refer to the pN and the J2 cases 
respectively. The pN values have been calculated considering 
an observer at L2 in the Sun-Earth system, and using the 
data reported in Murray and Dermott (1999) for the Solar 
System bodies. 
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well as their translational and rotational motion (table I). This implies 
that a suitable reduction of the GAIA observational data requires the 
development of a model of the celestial sphere more sophisticated than 
the one used for Hipparcos, the first global astrometry satellite launched 
by ESA in 1989. All mathematical and physical assumptions which 
characterize our model, as we shall specify in sections 2 and 3, are 
tailored to the GAIA concept. 

GAIA's main scientific objective is to shed definitive light on the 
origin, structure and kinematics of the Galaxy utilizing a complete 
census of the entire sky down to the visual magnitude V = 20, for a 
total of about one billion stars. This survey is achieved by exploiting 
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the same scanning strategy which was used with the successful astro- 
metric mission Hipparcos. In this way GAIA will measure the effects 
of many thousands of extra solar planets, and determine their orbit; 
thousands of brown dwarf and white dwarfs will also be identified. 
Many thousands of new minor bodies, inner Trojans, and even new 
Trans-Neptune objects, including Plutinos, may be discovered. The 
main physical properties of asteroids (10 5 — 10 6 new objects) will be 
investigated including masses, densities, sizes, shapes and taxonomic 
classes. In particular, asteroid mass is determined by measuring the 
tiny gravitational perturbations experienced in case of mutual close 
approaches (Bienayme and Turon, 2002). 

Moreover GAIA's fxarcsecond global astrometry allows one to test 
General Relativity (GR). Realistic end-to-end simulations in fact (Vec- 
chiato et al., 2003) show that GAIA could measure the PPN parameter 
7 to ~ 10 -7 (Id) after 5 years of continuous observations. The para- 
meter 7 measures the excess of curvature produced by mass-energy as 
compared to GR where its value is 1. As well known (Damour and 
Nordtvedt, 1993; Damour et al. 2002a, 2002b) a deviation from 7's GR 
value would deeply affect our understanding of fundamental physics. 

In what follows, greek indices run from to 3 and latin indices run 
from 1 to 3. 



2. Modeling the observations of GAIA 

The astrometric problem consists in the determination, from a pre- 
scribed set of observables, of the astrometric parameters of a star 
namely its coordinates, parallax and proper motion. While in classi- 
cal astrometry these quantities are well defined, in GR they must be 
interpreted consistently with the relativistic framework of the model. 

The first step of the modeling is to identify the gravitational sources 
and then fix the background geometry. Next we choose a suitable refer- 
ence system to label space-time points and describe the light trajectory, 
the stellar motion and the motion of the observer. Evidently the goal 
of our model is to write a formula which relates the observables to the 
astrometric parameters. To this purpose we had to take into account 
the way how the satellite GAIA will operate. 

Our first assumption is that the Solar System is the only source 
of gravity; in doing so we ignore, e.g., microlensing effects 1 which are 
perhaps the most important perturbations suffered by a light ray in its 

1 As far as the mission GAIA is concerned we expect ~ 25000 potential mi- 
crolensing events along the Galactic disk during the five year (Belokurov and Evans, 
2002). 
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way to us since they generate systematic errors in the data reduction. 
However, their expected number is very small as compared to the num- 
ber of stars surveyed by GAIA (> 10 9 ) hence we feel justified neglecting 
them at least at this stage of modeling. 

The second assumption is to consider the Solar System as a source of 
a weak gravitational field. This allows us to adopt a quasi- Minkowskian 
metric 

9a(S = Va0 + K(3 + 0(h 2 ) (1) 

where the h a p are perturbation coefficients such that \h a p\ <C 1; their 
spatial variations are at most of the order of \h a p\ while their time 
variations are at most of the order of (v/c) | h a g | where c is the velocity 
of light in vacuum and v is the typical velocity within the perturbing 
system. Notice that v represents the velocity needed to be gravita- 
tionally bounded in the Solar System (Damour, 1987), and its typical 
value for the energy balance (virial theorem) is ~ 10 kms -1 . Under 

these conditions we shall take h a p = J2 a wnere the sum is ex- 
tended to the bodies of the Solar System. The order of magnitude 
of each term can be expressed in term of powers of (v/c). Since 
the h terms are at least of order (v/c) 2 , the level of accuracy which is 
expected to be reached within the model is fixed by the order of (v/c) to 
which one likes to extend the calculation. Having these considerations 
in mind, we shall adopt the IAU recommended metric form for the 
barycentric reference frame inside the Solar System (IAU, 2000). In 
this case the lowest orders of magnitude of the metric coefficients are 
h 00 = 0((v/c) 2 ), h i:j = 0((v/c) 2 ) and h 0i = 0((v/cf). In the following 
we describe the construction of a model accurate to (v/c) 2 where in 
addition we approximate the Solar System to a static, non-rotating 
and non-expanding gravitating system and therefore we do not need 
to consider the spatial location of the individual gravitational sources 
at the corresponding retarded time. This assumption implies that in 
reconstructing the light trajectory of each photon emitted by a distant 
star, we keep the bodies of the Solar System fixed at the position they 
have with respect to its center of mass (CM) at the time of observation, 
say. Evidently each subsequent light ray will be reconstructed updating 
the positions of the bodies of the Solar System according to their actual 
motion. This does not mean that the bodies are at rest in the model, 
but that in our (v/c) 2 approximation their position can be considered 
fixed in a single integration. This model will serve as a touchstone for 
comparison once we extend our analysis to the order of (v/c) 3 . 

Let (r, x, y, z) be a coordinate system with respect to which the 
spacetime metric takes the form (1). Up to the order (v/c) 2 the space- 
time certainly admits a vorticity free congruence of curves which allows 
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for a spacetime foliation into three dimensional hypersurfaces with 
equation r =constant; we shall term them S(t) and identify the nor- 
mals with a unitary vector field u a = dx a /da with x a = r,x,y,z as 
a = 0, 1,2,3. In this case the spatial coordinates can be fixed within 
each slice up to spatial transformations only. We can then give the 
vector field u a the form u a = (— goo)" 1 '' 2 ^" which makes the vector 
unitary. The vector field so defined describes a physical observer with 
proper time a who Lie-transports the spatial coordinates. We also re- 
quire that the CM of the Solar System (figure 1) and the planets belong 
to this congruence. At this point of the modeling we have characterized 
the metric, the coordinate system and a physical observer u who is at 
rest with respect to the CM of the Solar System and therefore termed 
baricentric; we can now trace the light path. 




g= +h 

Figure 1. The astronomical set-up. A photon crosses the Solar System which is 
the only source of gravity. The vector field u describes physical observers who 
Lie-transports the spatial coordinates (including those of the planets at the {v/cj 2 
order). 
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Let T be a null geodesic with tangent vector field k a = d^ a /d\ 
(here A is a real parameter on T). A photon traveling from a distant 
star to the astrometric satellite within the Solar System would see the 
spacetime as a time development of the spatial hypersurfaces S(t). 
The light trajectory will end at the satellite at a coordinate time To 
and at a point with spatial coordinates x l Q . Because the spacetime is 
foliated into slices S(t) each representing a surface of simultaneity for 
the baricentric observer u namely her/his rest space, each point T(A) 
on the null geodesic can be mapped into the point on the slice S{tq) 
where the unique normal to S(t) at T(A) crosses S(tq), (see figure 2). 
In doing so the entire curve T is mapped into a curve Y in S(tq); this 
curve has a tangent vector £ a which at every point of Y is numerically 
equal to the components of the null vector k a in the rest space of u at 
the corresponding point of Y(\). 

Neglecting the time variations of the metric and imposing = 
the geodesic equation for the light ray leads, after some algebra, to the 
differential equations of Y. For the corresponding tangent field £ these 
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differential equations read: 

_ + e k ^-tdihoo) + ldih kj - -dkhijj tP - -d k h Q0 = o, (2) 

where £ is the unitary tangent vector field obtained reparametrizing the 
tangent field £, which is in general not unitary, and 5 hs is the Kronecker 
delta. 

To integrate equation (2) one needs to fix boundary conditions; 
these will be given by the components t of the vector field £ at the 
observation point, namely at r = tq and x % = x\. When the light ray 
is intercepted by the satellite then the angles between the direction of 
the incoming photon as seen by an observer comoving with the satellite 
and the three spatial directions of a triad adapted to that observer, will 
be the observables. These are essential to fix the required boundary 
conditions namely to express the components £\ in terms of known 
quantities alone. Our task now is to solve this problem. 

Let u' be the unitary vector field tangent to the satellite's world 
line and let {Aa} (where a = 1, 2, 3) be a space-like triad carried by the 
satellite. The expressions of the angle cosines e'a are given by a well 
known formula (see for instance de Felice and Clarke, 1990 or Brum- 
berg, 1991) containing the four-velocity of the satellite, the components 
of the triad which identifies the satellites's rest-frame, the metric coef- 
ficients at the observation time and of course the components i\ which 
are our unknowns. The most important terms entering that formula 
are the components of the satellite triad. As well known the latter 
identifies the rest-space of the observer comoving with the satellite. 
Such a triad is defined up to spatial rotations, namely transformations 
which leave the satellite's four-velocity u' unchanged. Each triad vector 
Aq, can be expressed in terms of coordinate components with respect 
to the coordinate basis {^ a } a =r,x,?/,z in the CM frame. If we denote the 
four-velocity of the satellite as 2 

u' = ( (d T + X s d x + Y s d y + Z S 8 Z ) (3) 

where ( is a function which makes the vector field u' unitary and 
X s , Y s , Z s are its x, y, z-components, a triad solution is given by 

Aj = Tid T + Xid x + Yidy 

Ag = T 2 d T + X 2 d x + Y 2 d y + Z 2 d z 

Ag = T 3 d T + Z 3 d z . 

where the T a , X a , Y a , Z a are the r, x, y, z-components of the correspond- 
ing Aa vector of the triad and are well defined functions of the com- 

2 The letter s is for "satellite" 
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ponents of the four- vector u' and of the metric coefficients, which are 
known quantities. As stated, we can relate the observed quantities e'a 
to the unknowns £^ by means of the above properties of the satellite 
frame. 

Although the general expressions of the boundary conditions that 
we have briefly mentioned will be published elsewhere, here we like to 
show them in the simple case of a satellite moving on a circular orbit 
around the Sun with coordinate angular velocity to. Recalling that the 
origin of the coordinate system is at the center of mass of the Solar 
System we obtain: 



nx 
Ho) 



e / i (x - x Q )^/g^V-n + (y - y^/fe^y/S + e'^-ffpo) 



ho) 



e'i(yo - y©) v / fev /Z n - (x - x )y^(w\/E + e'^-g 00 ) 



e' s yf-u. 



goo 



where xo,uq,zq are the coordinates of the satellite at the observation 
time, x Q ,y Q ,ZQ are those of the Sun, S = g X x{yo-yo) 2 +g y y{xo-x Q ) 2 , 
and II = c/oo + {g X x{yo ~ Ve) 2 + 9 yy (xo ~ x q ) 2 )uj 2 . 



3. Testing the model 

The extension of the many-body model to higher orders is expected to 
be complicated not only by the mathematical structure of the relativis- 
tic equations, but also by the numerical methods needed to implement 
it into a software code suitable for the reduction of the satellite data 
when they become available. Testing the model is of course essential; 
a basic step is to verify whether to the order (v/c) 2 and under the 
same conditions our model is consistent with the Schwarzschild non- 
perturbative model developed in de Felice et al. (1998; 2001). The 
tests that we have devised had the purpose to verify whether: i) the 
perturbative model restore full spherical symmetry when we consider 
the spherical Sun as the only source of gravity; ii) the light deflection 
caused by each individual body of the Solar System as evaluated in 
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our perturbative model agrees within the (v/c) 2 approximation with 
that expected in the Schwarzschild metric under the same observational 
conditions; hi) our model is able to reconstruct stellar distances. The 
tests as to the point i) were completely successful hence we shall only 
discuss those regarding points ii) and iii). 

An analytical formula for the light deflection as a function of the 
angular displacement ip of a star from the Sun, is given in Misner et al. 
(1973). Since that formula is itself approximated to the order of (v/c) 2 
and refers to a spherical Sun, we expect that its predictions coincide 
with ours, in case the Sun is the only source of gravity and spherical 
as well, at the {v/c) 2 order, that is at ~ 0.1 milliarcsec. Taking the 
same stars and computing the light deflection, the tests have shown 
a difference between the values deduced from our model and those 
obtained with the analytical formula of < 15 ^arcsec for almost limb 
grazing light rays. The difference becomes rapidly less than 0.1 //arcsec 
for tp > 5°. We are now performing a new series of numerical tests 
with our perturbative model to better investigate the reason of that 
discrepancy and the results will be published separately as soon as 
they are available (de Felice et al., 2003). 

In the exact Schwarzschild model we found a formula that expressed 
the components of the tangent vector k % to the null geodesic with 
respect of the spatial polar axes f, 6 and <p of a phase-locked tetrad 
associated to an observer moving on a circular orbit around the Sun. 
Given the position of a star and that of the observer (in Schwarzschild 
coordinates), we are able to deduce the cartesian components of the 
tangent to the null geodesic at the position of the observer by means 
of an almost completely analytical procedure, because the complete set 
of Schwarzschild coordinates of the star (r*, 8* ,(f>*) can be recovered by 
analytical integration. Practically, in this test we have considered the 
observer at two opposite positions on his orbit which we took as that 
of the Earth around the Sun (i.e. r\ = T2 = r®, 0\ = 02 = vr/2 and 
<t>2 = <t>i + Ti")- The stars are also lying on the orbital plane (9* = tt/2) 
so that the two observers are symmetrically placed about the Sun 
((f)* = (</>i + </>2)/2). The distance r* from the stars ranges from 1 pc to 
10 kpc. With this configuration it is easy to calculate the parallax of a 
star as p* = 1/r* rad (the distances being in AU), and so it is equally 
easy to determine the numerical accuracy of the position in terms of 
angles by simply subtracting the parallax p* of the Schwarzschild model 
and the approximate one reconstructed from the numerical integrations 
described above, namely p* a . The final tests show that the difference 
P* s ~ Pa i s ~ 18 ^arcsec for the entire distance range simulated. Here 
again we find a difference perhaps due to an insufficinet numerical 
accuracy and our present task is to improve on this situation. 
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